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Abstract 

The critical behaviour of the D = matrix model with potential perturbed by nonlocal 
term generating touchings between random surfaces is studied. It is found that the phase 
diagram of the model has many features of the phase diagram of discretized Polyakov's 
bosonic string with higher order curvature terms included. It contains the phase of smooth 
(Liouville) surfaces, the intermediate phase and the phase of branched polymers. The 
perturbation becomes irrelevant at the first phase and dominates at the third one. 
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1. Introduction 



After a recent progress, initiated by KPZ paper [0], in the investigation of conformal invariant 
field models with central charge c < 1 coupled to gravity, an understanding of physics beyond 
the "barrier" c = 1 is still an open question. In string theory, c > 1, KPZ approach leads 
to complex anomalous dimensions of fields indicating that an instability is developed in the 
theory. 

The same phenomena was noticed independently under computer simulations of discretized 
Polyakov's bosonic string. It was found 0] that near the point with c = 1 transition to the phase 
of crumpled surfaces occurs. In this phase the world sheet of string has a form of "hedgehogs" 
and an underlying effective theory cannot be a string theory. One of the well known ways to cure 
this instability is to perturb the string action by adding higher order intrinsic curvature terms, 
like / (Pa^/gR"^. It turned out that the discretized perturbed model has four different phases 
depending on the values of central charge and constant coupled to the intrinsic curvature: phases 
of crumpled surfaces, branched polymers, intermediate phase and phase of smooth (Liouville) 
surfaces. It is in the last two phases where one may hope to construct an effective string theory 
with the central charge — oo < c < 4. For —oo < c < 1 the model has a remarkable universality 
property: additional curvature terms become irrelevant in the continuum limit and the critical 
behaviour is governed only by the central charge, in accordance with KPZ relations leading to 
negative values of string susceptibility exponent ■jstr < 0. At the same time in the intermediate 
phase, for 1 < c < 4, the curvature terms are responsible for nontrivial critical behaviour of the 
model with positive string susceptibility exponent 

< -fstr < 0.5 (1.1) 

An understanding of this phase in the continuum theory is still lacking although possible mech- 
anism of its appearance was advocated in It was proposed [^, || to use the matrix models 
to justify the critical behaviour (|0|). 

Matrix models allow us to reproduce KPZ results for conformal matter with c < 1 0. 
In order to get the critical behaviour ( |1.1| ) it seems natural to perturb the potential of the 
matrix models by terms which include the effects of higher order curvature terms. Possible 
form of such terms was discussed in [|]. The simplest one was found to be (TrM^)^ and 
possible matrix model, D = hermitian one-matrix model with the potential Tr (M^ + gM^) 
perturbed by this term, was considered in 0]. The appearance of additional phase with positive 
exponent -^str = 1/3 was found and the existence of a series of analogous phases with '~^str = l/'"-, 
(n = 3, 4, . . .) was conjectured. Although these points lie in the region ( [l.lD their interpretation 
as critical points of c > 1 matter coupled to gravity remains unclear. 

In the present paper we investigate the influence of the perturbation (TrM^)^ on the critical 
behaviour of = hermitian one-matrix model with an arbitrary polynomial potential. We 
identify the critical points of the model and find that the resulting phase diagram has many 
features of the phase diagram of the discretized Polyakov's string with higher order curvature 
terms included Under increasing of the constant coupled to perturbation the system passes 
through the phase of smooth surfaces, intermediate phase and the phase of branched polymers. 
It turns out that the perturbation is irrelevant at the first phase and dominates at the last 
one. We show that in the continuum limit the intermediate phase with positive exponent 
^str = + 1) {K > 2) appears only as perturbation of (2, 2K — 1) minimal conformal model 

coupled to gravity. 
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2. Potential of perturbed model 



The partition function oi D = hermitian matrix model is defined as 



J dM exp{-aNTTV{M)) (2.1) 



wliere a is tlie bare cosmological constant, V{M) is the potential of the model and integration 
is performed over hermitian N x N matrices M. For large N and a close to the critical value 
ctcr, at which the perturbation series for Z{a) become divergent, the string susceptibility has a 
behaviour 

with nontrivial critical exponent 7^4^- It is well known 0, that m. D = Q hermitian one-matrix 
model the string susceptibility exponent can have two different kinds of values associated with 
multicritical Kazakov models and Penner model [§, |^ . In the K—i\i multicritical model the 
potential has the following form 

V{M) ^ UMM) = (^y=^^^M- (2.3) 

and the exponent •jstr takes negative values 

Istr = (2.4) 

This model was identified in the continuum limit as (2, 2K — 1) minimal conformal model 
coupled to 2D gravity [0]. In the Penner model with nonpolynomial potential 

y(M) = log(l - M) + M (2.5) 

the critical exponent is equal to 

7str = (2.6) 

and in the continuum limit model possesses many properties of c = 1 quantum gravity |^, |^. 
In both cases (|2.4|) and ( p.6|) the exponents 'jstr are in agreement with KPZ predictions. 

Let us perturb the potentials ( p.3|) and (|2.5|) to include higher order curvature terms. The 



simplest way to do it is to add the following "nonlocal" interaction term (TrM^)^ to the 
potential [|]. To this end one chooses the potential of perturbed matrix model as 

Trl^(M) = Tryo(M) - ^{Ti (2.7) 

where Vo(M) is the generalized Penner potential^ 

TrK)(M) = J dxp{x)TT log(x - M) (2.8) 

^The matrix model with potential Vo{M) was considered in pl| where its nontrivial critical behaviour was 
found provided that one identifies cosmological constant with the dimension of the matrix AI. 
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and p{x) is the normalized function with well defined moments 



J dx p{x) = 1 j dxx ^p{x) = finite 



having a meaning of the density of external field. The reason for this particular choice of the 
potential is twofold. First of all, expanding potential Vo(M) in powers of M it is possible to get 
an arbitrary polynomial potential by tunning of the moments / dxx~^p{x). Secondly, it turns 
out that the presence of an external field p{x) simplifies the solution of the model. 

To understand the meaning of the perturbation (TrM^)^ one considers the form of random 
surfaces generated by the potential (|2.7D. The first term, Vo(M), being expanded in powers of 
M leads in a standard way to a sum over discretized (in general unconnected) closed surfaces 
of an arbitrary genus. The last term in (|2.7| ) opens a possibility for these surfaces to touch each 
other. Then, unconnected surfaces with touchings included become connected and have a form 
of "soap bubbles" [Q. Note, that the touchings may change the genus of the surface. 

In what follows we consider the model in the leading N ^ oo (spherical) approximation 
which allows us to identify the possible critical points. In the spherical approximation any 
two planar surfaces may touch each other only at once and cannot form closed chains of soap 
bubbles. A power counting indicates that the partition function ( p.l|) has the following 
decomposition for large 



(a) = ^ Z(^, r)a- V + C(iV°) (2-9) 

S,T 



where S is the total area of the discretized planar surface, T is the number of touchings this 
surface has and Z{S, T) is the partition function for fixed S and T. This expression allows us 
to identify the bare cosmological constant as a = e^-^ > 0. In the continuum limit when a 
approaches the critical value c^cr the surfaces with area — > oo dominate in the sum (|2.9|) with 
the following asymptotics Y.t Z{S,T)g^ ~ a^j,{g)S'^'"'''^^^~^, leading to the critical behaviour of 
the partition function Z{a) ~ (a — acr{g)Y~'^^''''^^^ ■ Here, the exponent 'jstrid) and the critical 
value of cosmological constant are functions of touching coupling constant. Formally, acr{g) is 
given by series in g and nonanalyticity of this function for some g would mean that the critical 
behaviour of the model is determined by the surfaces with infinite number of touchings. One 
expects that in the limit acr{g) ^ g when it is more favorable for the critical surfaces to form 
a touching instead of increase area the critical surfaces are degenerated into infinite chains of 
"soap bubbles" and the model turns into the phase of branched polymers |^ . 

To find the phase diagram of the model in (a, g) plane we calculate the string susceptibility. 
Using the definition (^.2[) and an explicit form of the potential one gets 




X = dxpix) (^Tr log(x - M)) - f ((^TrM^ 

where for an arbitrary operator {A{M)) = e"^'-"-* / dM A{M) exp{—aNTTV{M)) denotes the 
connected correlator. Here, the second term is responsible for the touchings of random surfaces. 
It can be simplified for large using the factorization property of the vacuum averaged product 



of invariant operators |12 



((TrM2)2) = ((TrM'))2 + C(ArO) ^2.10) 
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This property implies that in the spherical approximation the touchings of the surfaces are 
measured effectively with the following constant 

-g = -ga(^^i:iM^'^ (2.11) 

called effective touching coupling constant. Then, for the string susceptibility we get 



X 



J dx p{x) ^^Tr log(x — M) 



da \J \N I 2ga ) 

The explicit form of this expression suggests us to connect x with one loop correlator 



z-M 



using a simple identity 

/ 1 



dz {W{z)) = l^—Ti log(x - M)^ - log A, 



as A ^ oo 



Here, the divergence of the integral as A — oo is a consequence of the asymptotics of loop 
correlator 

{W{z))-^l/z, as^^cx). (2.12) 
Finally, the expression for string susceptibility is given by 

^ = -^[jd^ P(^) jl dz {W{z)) + log A - . (2.13) 

One completes this expression by the relation between effective touching coupling constant 
( |2.11| ) and one loop correlator 

g = -gaj^^^z^W{z)) (2.14) 

where integration contour C encircles singularities of {W{z)). Equations ( ^.13[ ) and ( P^.14| ) allow 
us to find the string susceptibility starting from one loop correlators. 



3. Perturbed loop equation 



To calculate loop correlators one uses a powerful method of loop equations 0, ^ [T3|, [T^. In 



D = Q hermitian matrix model with potential U{M) the solution of loop equations gives us 
(one-cut) loop correlator in the limit — > cxd as ||I4[ 



2Jc 2m z - u> ^(^-A+)(^-A_) 

where boundaries X± of the cut are fixed by the asymptotics ( |2.12| ). Since the potential of 
the model (p.7|) contains an additional nonlocal term the modification of the loop equation is 
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needed. The perturbed loop equation we get in a standard way |T3[ performing the change of 
the integration variables M — ^ M + in (|2.1| ) with infinitesimal e 



{W{z)W{z)) = a 



duj 1 
c 2TTi z — uo 



{V^{uj)-J-^TtM')W{u] 



Using relation analogous to ( p^.l4D it is easy to express Tr as an integral of W{z') and then to 
get a closed equation for loop correlators. In the limit iV — > oo we apply factorization property 
( 2.101 ) to transform the loop equation into 



{W{z))^ = a 



duj U\uj) 
c 2m z — ui 



where the notation was introduced 



U\z) = V^{z)-zga{-^TiM' 



dxp{x) 



X — z 



+ zg 



(3.2) 



(3.3) 



The solution of the equation ( p.2| ) is given by ( ^.1| ) provided that one considers U'{z) as a 
derivative of some potential. Let us choose for simplicity the density p{x) to be an even 
function 

p{-x) = p{x). 

Then, the potential V{z) is also even and, as a consequence, the cut of the loop correlator is 
symmetric with A+ = — A_. After substitution of ( |3.3| ) into ( |0| ) the loop correlator is equal 
to 



{W{z)) = - 



a f dx p{x) 



1 



-)j + ^g{z - Vz^ + c) 



(3.4) 



z - X \ \fx^ c 

and the additional condition (|2.12|) leads to the following equation for c = A+A_ = — A^ < 



-g = 2/(c) - — 
ac 



Here, the notation was introduced for the function 

dx p{x) 



\/a?M-~c {x + \/ x^ + c) 



- [ dxpix) I 1 , 1 



(3.5) 



(3.6) 



depending on the moments of density p{x). For latter convenience one defines the following 
functions 



^(c) 



dx p{x) 



(X + V^^T^)2' 

They are related to each other by simple equations 



0(c) = c 



dx p{x) 



yxM-c (x + a/x^ + c)' 



or equivalently 



/'(c)+0'(c) = — 0(c), 
c 



dtrfit) 



/(c) + 0(c) = 2^(c) 



c^ Jo 



^(c) 



C^ Jo 



dttf{t) 



(3.7) 
(3.8) 



Note, that loop correlator ( p.4| ) depends on the effective touching coupling constant which can 
be found in its turn after substitution of {W{z)) into (|2.14|) as 

lQ + g[acyj + c (x + Va;^ + c)2 IQ + g[acY 



Being combined together equations ( p.4| ), (|375|) and ( |3.9| ) define one loop correlator in the 
spherical approximation. 



4. String susceptibility 



To evaluate the string susceptibility we substitute one loop correlator 
After integration one gets in the limit A ^ oo 



into relation ( |2.13D . 



/ dxp{x) £dz{W{z))+ log A =^JdxJ dyp{x)p{y)\og f^^y±^±^/^±W¥T^ 



+ (1 — a) J dx p{x) log 



X 



^ ~^^9 J dx p{x)x{x — y/x'^ + c) + ^gc (4.1) 



As follows from (|2.13|) , the string susceptibility contains a full derivative of this expression 
with respect to a. Its calculation is not trivial because both c and g depend on a. However, 



expression ( [4.1|) has a remarkable property: its partial derivative with respect to c equals zero 
as only one takes into account condition ( p.5| ). Hence, differentiating ( |4.1| ) we may ignore the 
dependence of c on a and after some transformations 



X = I dx I dyp{x)p{y)\o^ 



ga da \ 8g 



( 



\ 



i + ^ + yrr 



1 + 



g{acfi){c) 



J 2ga'^ \~ ' Ag 
Here, g is given by ( p.9|) and the equation for c we find from ( ^75|) and (|3.9| ) as 

2{lQ + g{acy) 



(4.2) 



ac 



16f{c)+g{acy^{c) 



(4.3) 



Solving the system of equations (^3), ( |4.3|) and ( p.9|) one can evaluate the string susceptibility 
X for an arbitrary even density p{x). The properties of x and possible critical behaviour of the 
model are considered in the next section. 



5. Analysis of critical behaviour 

It is convenient to start the analysis of critical behaviour by considering the special case g = 
when the potential ( p.7| ) reduces to the potential ( |2.8| ) of so called "reduced" hermitian matrix 
model with even density p{x). 
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5.1. Unperturbed critical behaviour 

For g = the equations (|4.3| ) and are replaced by 



- = c/(c) g = (5.1) 

a 



The string susceptibihty is given by the first term in (O) and depends on a only through the 
function c = c{a). Therefore, there appear the following two sources of nonanalyticity of the 
string susceptibility x = : singularities in the dependence of x on c (Penner's point) 

or in the dependence of c on a (Kazakov's multicritical points). For instance, for even density 
p{x) having a form 

p{x) = l{5{x-ri) + S{x + ri)) (5.2) 
the string susceptibility (|4.2| ) has a logarithmic singularity at Ccr = —rf. The corresponding 



critical value of the cosmological constant can be found from (^.1|), (^.2|) and (|3.6|) to be a, 



^ ^ = and near this point a = —2^1 + Hence, the string susceptibility (|4.2|) has a 
logarithmic behaviour ^ 

X = -|loga (5.3) 

as a — > 0. Such a behaviour of x is one of the features of c = 1 quantum gravity. Note, that 
despite of the original Penner model 0, ^ the critical point a = 1 did not appear here because 
the underlying mechanism of criticality based on the shrinking into point of the boundaries 
of the cut of one loop correlator (A+ — > A_ 7^ 0) is suppressed for even potential Vb(M) 
leading to the symmetric cut with A+ = — A_. Penner critical point originates from the singular 
dependence of x on c. Otherwise, for x having a well defined expansion in powers of c, the 
nonanalyticity of the function c = c{a) gives rise to criticality. In particular, by choosing the 



function /(c) defined in (|3.6|) as 



cf{c) = 2{l-[l + ^y)=^ (5.4) 



one gets c — Ccr = 4(a — OcrY^^ for acr = 1 and Ccr = —4. Then, the differentiation of ([4.2| ) 
yields 

dx 
dc 



3=0 



and near the critical point (a —>■ acr) the string susceptibility is equal to 

X = -\{C- Ccr) = -(a - C^cr)^^^ 



leading to the i^— th multicritical point with the exponent ( p.4| ). By comparing the explicit 
form ( ^.4|) of the function /(c) with the definition (|3.6| ), we find the moments of the density 
p{x) and restore the corresponding potential Vo(M) to be equal to 



5.2. Perturbed critical behaviour 

For g ^ when touchings of the surfaces become to play the role the string susceptibility 
( |4.2| ) acquires the additional term depending on the effective touching coupling constant g and 
its derivative. As a result, there appears a new mechanism of criticality additional to that 
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considered before for g = 0. This mechanism explores the singularities in the dependence of ^ 
on a and leads, as it will be shown below, to the appearance of new critical points with positive 
exponents ■jstr- 

We begin the analysis of the system of equations ([4.2|), and (|3.9|) by calculating the 



derivative The zeros of this derivative would mean the nonanalyticity of the function c 
c{a) which is one of the sources of criticality. Differentiation of the both sides of ( [4 .31) gives 

da oP'c 



dc 2(16-c/(ac)2(l-ac0(c))) 



(r(c) + — )(16-,M^) (5.5) 



In the special case g = this expression with (|5.4| ) reduces to 



da 
dc 



k c\^-i 

1 + - , as c Ccr = —4 



n 4 V 4 

9=0 



and the derivative has K—th. order zero. One finds that for g ^ the equation ^ = has two 
different branches of solutions: 

\2 



and 



g{acrCcr) = 16 (5.6) 

/'(ce.) + = (5.7) 

It is interesting to note that the first equation does not depend on the explicit form of the 
density p(x), whereas the second one does not depend on the touching coupling constant. 
Hence, the appearance of the (|5.6| ) is a general property of a matrix model perturbed by 
touchings interaction term. It is this solution that leads to the phase of branched polymers for 
g > go when only touching term is responsible for the formation of the critical surfaces and 
an explicit form of the potential Vq{M) turns out to be unessential. At the same time, only 
solution ( ^.7|) survives in the limit g ^ and it corresponds to the phase of smooth surfaces. 

As we saw before, in the unperturbed matrix model with potential Vo(M) the existence of 
solutions of the equation ^ = leads to negative critical exponents jstr for which the string 
susceptibility ( |2.2D is regular as a — > acr- In our case, (? 7^ 0, the expression ( |4.2| ) for x contains 
the derivative 




which tends to infinity for the solutions (|5.6|) and ( ^.71) of the equation ^ = provided that 
^ 7^ 0. An important property of the model is that the derivative of the effective touching 
coupling constant ^ is finite for the solution ( |5.7| ) and infinite for ( |5.6| ). Indeed, using the 
expression ( ^Tsf) for g we obtain 




and after substitution of the expression ( [5. 5] ) for ^ 

dg 4:gac^(j){c) 
da 16 — g{acy 



(5.8) 
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This relation enables us to identify the peculiarities of the derivative Firstly, equation ( |5l8| ) 
is singular for (|5.6|) and regular for dSTT]). Secondly, ^ depends on the function 0(c) and it 
would have additional singularities if this function is singular for some c as it occurs for c = —rf 
in the Penner critical point with the density ( p.2| ). 

We forget for a moment about a second possibility and return to its discussion at the end of 
the paper. In that case, the singularities of the derivative ^ and zeros of ^ become correlated 
and there are only two branches of criticality: (|5.6|) and ( |5.7|) . For the solutions ( |5.7| ) we have 
^ = and ^ = finite. It means that the string susceptibility ( [4. 21) is finite for a — > and 
the critical exponent is negative '-jstr < 0. For the solutions (|5.6| ) the relations ^ = 0, ^ — cxd 
lead to the infinite string susceptibility as a — *• acr- The critical exponent becomes positive and 
we expect to get the value 'jstr = 1/2, corresponding to the phase of branched polymers. It 
turned out that these two branches may cross each other given rise to intermediate phase with 
positive critical exponent. To conclude this section we stress that it is a nonanalyticity in the 
dependence of effective touching coupling constant on the cosmological constant which leads to 
the appearance of phases with positive critical exponent '-fstr- 



6. The properties of critical points 



In this section we consider the properties of the solutions ( ^.6|) and (|5.7|) . The equation (|5l7| ) 
relates the critical values a = acr and c = Ccr for which ^ = and 



a — acr ~ (c — c, 



2 

cr J • 



To value of Ccr one gets after substitution of (|5^ ) into (|4.3|) as a solution of the following 
equation 

4:{Ccrf'{Ccr)f{Ccrf{Ccr)y = giCcrfpiCcr))' ■ (6.1) 



The effective touching coupling constant at the critical point is found analogously from (|3.9|) , 
d^) and (§3) as 

g{acr) = 2{Ccrf{Ccr))' (6.2) 

and its derivative (|5.8|) is given by 



d9_, . _ 2gcl^4>{Ccr)f'{Ccr) , . 

da^""-^ 4{Ccrf'{CcrW-g ^ ' 

These expressions are finite for all values of Ccr except of the point where {ccrf'{ccr)Y = \g 
and the equations ( ^.6| ) and ( |5.7| ) coincide. With all these expressions taking into account we 
obtain from (|4.2|) the string susceptibility as x ~ c — Ccr ~ (a — ctcr)^^^- In general, by tunning 
of the density p(x), it is possible to satisfy ^ ~ (c — Ccr)^^^ in order to reach the multicritical 
points with x ~ (a — acrY^^ and negative exponent 7^4^ = —^/K. 

Considering the second branch ( p.6| ) we note that it exists only for the positive touching 
coupling constants (? > 0. Since the equation (|5.6| ) does not depend on the density p(x) it is 
impossible in contrast with the previous case to set higher derivatives n > 1 to zero at the 
critical point and the only possibility remains 

a - acr ^ {C - Ccrf 
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The critical values Ccr and are found from ( [5.6| ) and ( [4.3| ) as solutions of the system 

4 



/(Ccr) + 0(Ccr) = "V^ 



The effective touching coupling constant (|3.9| ) has a form 

(7(acr) = /(Ccr) - 0(Ccr) = 2(CcrV'(Ccr))' 

and is finite at the critical point, whereas its derivative ( |5.8| ) has a pole as c 



c - 



(6.4) 



(6.5) 



(6.6) 



Hence, being substituted into (^7^) , this derivative defines the asymptotics of the string sus- 
ceptibility near the critical point as 



X 



dg 



Ail){ccr) da 



(a — a,- 



-1/2 



(6.7) 



The critical exponent jstr = 1/2 corresponds to the phase of branched polymers. 

The solutions of the equations (|0| ) and (|6^ ) define two curves of the possible critical 
behaviour in the (a, g) plane. As was shown in @] in the simplest case of potential Vo(M), 
these curves are intersected at some point leading to appearance of new intermediate phase. 
The condition for the existence of such point can be found by combining equations (|5.7|) and 
into the system 



/(Ccr) + 0(Ccr) = -^9, 



Ccrf (Ccr) 2V^' 



(6. 



The critical point defined by these equations belongs to both branches of solutions and has 
the properties of both. In particular, relations ( |6.3| ) and ( |6.6| ) coincide for the solution of (|6.8| ) 
indicating that the derivative ^ has a pole at the critical point and, as a consequence, the string 
susceptibility has asymptotics x 
equal to zero for 



Moreover, the both factors in the expression ( |5.5| ) for 



da 
dc 



As only ^ ~ c — Ccr on the branch ( |5.6| ) and ^ ~ (c — Ccr)^ ^ 



dc '^^ N — —J' dc 

after tunning of parameters of model on the branch ( [5.7[ ), then at the intersection point one 



gets 



do. 



dc 

X ~ (a — ftcr)^"*^ with positive exponent '-jstr = l/(-^ + !)• 

We conclude that the existence of solutions of the system ( |6.8|) would mean that between 
the phase of multicritical points with 7^*,- = —1/ K and phase of branched polymers with '-^str = 
1/2 there appears an intermediate phase of "hybrid" multicritical points with the exponent 

^str = l/iK+l). 



(c Ccr) 



K 



or a 



ar 



[C-Cr 



)^^^. Hence, the string susceptibility has a behaviour 



7. "Hybrid" multicritical points 

The necessary condition to get the hybrid multicritical point is to satisfy the relation a — ac 
(c — Ccr)^ near the critical point defined by equations (|5^ ) and (|6.1|) , or equivalently. 



da 
dc 



d'^a 



d' 
dc^- 



a 







(7.1) 
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At the vicinity of the critical point ( |5.7] ) the derivative ( |575| ) is given by 



dc ac 



and condition ( [7.1|) can be satisfied under an appropriate choice of function /(c). Namely, 
differentiating the both sides of this relation with respect to c and using ( 171] ) 



we obtain set equations for the function /(c 



/"^ce.) + ^C^' = 0' (n = l,...,ir-l) 



The remaining higher derivatives are not fixed and the function /(c) is defined up to 0{{c 
Ccr)^) terms. The general form of this function is 



9 / / r\^\ 



/(c)=^ + — 1- 1-- +O((c-co)") (7.2) 
where arbitrary constants A and ao are defined as 

2 

A = /(co) + co/'(co) ao = ~^|j7^ C^-^) 
and the critical values of c and a are given by 

Ccr{g) = Co, acr{g) = "o (7.4) 

The expression ( |7.2|) has a well defined expansion in powers of c with the coefficients related 
to the moments of the density p(x). Using identities ( ^.8] ) one can evaluate the functions 0(c) 
and ip{c), substitute them into equation ( |6.1| ) to find 

g = 4{K + l)- ^ 



ttoco 



This equation implies that the constants A, and Co, as well as the critical values of Ccr 
and acr are functions of touching coupling constant g. After their substitution into ( |7.2| ) the 
function /(c) also becomes g dependent. However, according to the definition ( p.7|) of the 



model the density p{x) and function /(c) do not depend on g. The identity = leads to 
additional restrictions on the possible form of /(c). After substitution of ( [7.2|) this equation has 
the following two solutions 

dcojg) _ dao{g) _ dA{g) _ ^ 
dg dg dg 

for an arbitrary K > 2 and the second solution exists only for the special value of K 

dA{g) 4 daoig) 2ao{g) 



dco ao(5')co(5')' dco Co{g) 



for K = 2 (7.6) 
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For the last solution /(c) is a linear function of c and the corresponding potential has a form 
Vo(M) ~ + aM'^ with some a. It is this model that was considered in @]. Notice, that 
in order to get multicritical points with K > 2 one has to consider only first solution ( |7.5| ). 
Relation ( [7 .51) implies that on the branch ( ^.71) of criticality the critical values Ccr and acr do 
not depend on the touching coupling constant g. On the other hand, the value of Ccr one gets 
solving equation ( |6.1|) which depends on g. The only way to satisfy the both conditions is to 
set 

(c/(c))'L=,„ = Co/'(co) + /(co) = 

or to put y4 = in expression ( [7. 2| ) for /(c). The resulting function /(c) coincides (for ao = 1 
and Co = —4) with an analogous function ( p.4|) leading to Kazakov's multicritical points for 
g = 0. However, for g ^ equation ( |6.1|) leads to an additional constraint on the density p{x) 



(C0(c))'|^^^^ = Co0'(Co) + 0(Co) = 

or using (|3.7]) 



/(Co) = 0(co) /'(Co) = 0'(co) (7.7) 

It can be easily verified that these conditions are not fulfilled if one chooses for /(c) the ex- 
pression (|7.2| ) (with A = 0) without 0{{c — Cq)^) terms. The "minimal" expression consistent 

with the condition ( |7.7| ) one gets by adding to /(c) the following term: ^^^^^^ (l — ^) and 
after simple calculations 



c/(c) 


2 








2 


c</)(c) 










2 








(7i 



Using the definition ( |3.6| ) of the function /(c) we calculate the moments of the density and then 
find the potential Vo(M) as 

yo(M) = -— Uk + 1)U2k{M') -{K + 2)U2(k+i){M')) (7.9) 

where M' = -^^^M and U2k{M) is the potential ( |2.3| ) of the K—th multicritical Kazakov's 
model. It is this potential which leads to the appearance of the hybrid multicritical points. 

Let us consider the critical behaviour of the model with the potential ( [7.9| ) following the 
general analysis of the previous section. 

Smooth surfaces: g < go 

For g = the potential ( p.7|) is given by ( |7.9| ) and the model has a critical behaviour corre- 
sponding to the i^— th multicritical point. The values of the critical parameters are given by 
(|7.4| ) and they lie on the branch ( |5.7| ) of solutions. Moreover, for g < this branch is the unique. 



After increasing of g > the model follows the same branch as for = until the touching 
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constant reaches the critical value at which the both solutions ( p.6| ) and ( ^.Tp coincide. The 
value go is found from the system (|678|) and ( [7.8|) as 



16 

^° (aoco)^ 

Thus, for g < go the critical behaviour of the model is described by the equations ( ^.71) , ( |0| ) 
and (T^). This phase has the following remarkable properties. It follows from ( |7.4| ) and ( |7.5| ) 
that the critical values ctcr and Ccr don't depend on g and coincide with analogous values (|7^ ) 
for the i^— th multicritical model without touchings. Substituting ( [7.4| ) and ( |7.8| ) into ( |6.2| ) we 
find that the effective touching coupling constant vanishes at the critical point 

giacr) = 0, for g < go (7.10) 



To get the behaviour of a and g near the critical point (|7.4| ) one evaluates their derivatives 
(|5.5| ) and and gets 



l-^)=(A-+l)*^(l-£)"+0((c-c„) 



A'+l 



and 

^(«) = -^^^(« - ao) + O ((a - ao)"^) (7.11) 

^ go- g ^ ^ 

Note, that (7 is not analytical in a near the critical point ctcr = «o and analytical in (yf for g < go- 
It means that at the vicinity of the critical point the effective touching constant g is defined by 
the random surfaces with infinite area and only finite number of touchings. One uses ( [4.2| ) and 



2 



( |7. 1 1| ) to calculate the derivative = (f^j Hence, the string susceptibility 

(|4.2| ) in this phase is given by 



X ~ (c - Co) ~ (a - aoY^^ 
indicating that for g < go the model looks like ordinary K—th multicritical model. 

Intermediate phase: g = go- 

For this value of g the solutions ( |5.6| ) and ( |5.7| ) coincide giving rise to new critical point with 

16 

Ccrigo) = Co, acr(fl'o) = «o 5' = fi'o = ^^^^^y C^-l^) 

To find the scaling of parameters near the critical point one substitutes ( |7.(j| ) and ( |7.12| ) into 
(|5.5| ) and obtains 



dc Co \ Co 

and integrating 

« \ _ ^ C \ 
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B 


as 




K 


a ^ 


y K + 1 







g(^a) = ^{K +l)K-T^ (l )'"■" (7.13) 

aoCo V aoJ 

Although it vanishes at the critical point ( |7.12| ), the derivative ^ ^ oo tends to infinity. The 
asymptotics of the string susceptibility (|4.2|) in this phase is determined by the singularity of 



the derivative ^ 

da 

X = -^«oC^^(co)£~(«-«o)-^/(^+^^ 
with positive critical exponent jstr = l/(-^ + !)• 
Branched polymers: go < g < g-max 

For large values of the touching constant the model passes to the phase of branched polymers 
described by the relations (|6.4]) and (|7.8|) . The equation (|6^ ) for critical point looks like 



Co / Oicrig) 

and after it substitution into (16.41) and 



1 / 1 

go\ ' acr{g) / , /, ao \ 



\g J ao \ V oicrig), 

and 



(7.14) 



K 



AjK + l) ^ ao 
g{acr) = 1 ^ 

These relations give us the dependence of the critical values of cosmological constant ctcr and 
effective touching coupling constant g{acr) on g as shown in fig. 1. We notice that for g > go the 
critical point starts to depend on the touching constant. In the phase of branched polymers the 
critical value ttcr is increasing function of touching constant and in the limit when g approaches 
a maximum value the cosmological constant tends to infinity 



g gmax = go{K + l) , acr{g) oo, g{acr) 

aoco 

In an opposite limit, g go, the dependence of acr and g{acr) on g is found from ( [7.14| ) as 
(ocM \ ^-K-i(g .Y'^^ X K+l ( g Y 

These relations imply that under transition of the system through the critical point with g = go 
and acr = «o the derivatives and ^^y_f undergo a jump. 

Let us consider the behaviour of parameters near the line = c^crig) of criticality. Since 



da 
dc 



and 



d^c 
dc- 



7^ then a — Ocr ~ (c — Qr)^. As follows from ( |6l6D , ^ 
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Vg Ccr4'iCcr) ^ 1 
2 C—Ccr ^ OL — OLcT 

relation gives 



and the derivative of g tends to infinity as a — > acr- Integration of this 



g{0L) - g{acr) ~ V" - "cr- 

As was mentioned before, the string susceptibihty in this phase is equal to 
exponent jstr = 1/2 corresponding to the branched polymers. 
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(7.15) 
with the 



90 



ao 



a 




(a) 



(b) 



Fig. 1: The phase diagram of the matrix model with perturbed polynomial 
potential. It contains the phase of smooth surfaces (g < go), intermediate 
phase {g = go) and phase of branched polymers {g > go)- (a): The dependence 
of the critical value of cosmological constant on the touching constant; (b): 
the effective touching constant along the line of criticality. 

We found that the model has three phases with the string susceptibility exponents •jstr = —l/K 
for g < go, jstr = ^/{K + 1) for g = go and jstr = 1/2 for g > go- The position of the critical 
points in shown in fig. 1(a). Let us consider the properties of the random surfaces contributed 
at the partition function in different phases. The critical surfaces have an infinite area and may 
touch each other. The average number of touchings (T(a)) = we obtain in the large 

hmit using the definitions (|2l|), (|2lO| ) and (PD as 



(T(«)) 




TrM^ 



a 
T 



\n 



^9 



2 + 0{N-^) 



The dependence of the effective touching constant on g is shown in Fig. 1(b). Using it we find 
the average number of touchings as 



(T(acr.)) = ioi g < go. 



(T(«cr)) > ioi g <go. 



Thus, at the phase of multicritical Kazakov's points the critical surfaces have not touchings 
and the perturbation becomes irrelevant. It is interesting to compare this property with an 
analogous one in the model corresponding to the solution ( |7.6| ) for K = 2. By choosing /(c) 
as a linear function of c and proceeding through relations (|b.lD, ( |5.7| ) and ( |0| ) one gets that 
in contrast with the previous case both the effective touching constant g and touching number 
{T{acr)) differ from zero in the phase of smooth surfaces. 

The behaviour of (T(a)) near the critical line of fig. 1(a) can be found using the asymptotics 
(|7.11| ), (|7.13| ) and (|7.15|) of the effective touching constant: 



9<gc 



{a -acr)' 



{T{a))\ 



a=9c 



(a-acr)^+i 



{na))\ 
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These relations imply that under transition through the point g = go the second derivative of 
the partition function 772 '^dadg' ~ ^ undergoes an infinite jump. 

There are two ways for critical random surfaces to increase the area: to form an additional 
"soap bubble" or to increase of the area of the existing ones. At the phases of smooth surfaces 
and branched polymers only one of these mechanisms works. At the intermediate phase one 
expects that both mechanisms are responsible for the formation of critical surfaces. 



8. Conclusion 

Let us summarize the properties of matrix models perturbed by "higher order curvature term" 
(TrM^)^ generating the touchings between the random surfaces. We found that the matrix 
model with perturbed polynomial potential has three phases: smooth surfaces, branched poly- 
mers and intermediate phase. At the intermediate phase which appears for special value go 
of touching constant g the string susceptibility exponent has positive value •jstr = ^/{K + 1) 
for some K > 2. For g < go the model turns into the phase of smooth surfaces corresponding 
to K—th multicritical point. We obtained that the critical surfaces haven't touchings. The 
perturbation becomes irrelevant in this phase and the model is equivalent to (2, 2K — 1) min- 
imal conformal model coupled to gravity. For g > g^ the critical behaviour is dominated by 
the perturbation. The random surfaces are degenerated into branched polymers and the string 
susceptibility exponent approaches the maximum value 'jstr = 1/2. 

In conclusion, we briefly discuss the critical behaviour of the perturbed matrix model with 
Penner like potential defined in ( |2.8|) and ( p.2|) . As was shown in sect. 5.1, the singularities 
of the function /(c) at c = —1]'^ play the central role in an analysis of unperturbed model. 
The same is hold also for g 0. We have two sources of criticality in ( |4.2| ): logarithmic 
peculiarities of x for c = —rj'^ and singularities of the derivatives ^ which define two different 
branches of possible critical points. At the first case, the singularities of the functions /(c) 
and 0(c) lead in ( [4.3|) to = for an arbitrary touching constant g including g = and 
string susceptibility has logarithmic behaviour (|5.^). At the second case, the derivative ( |^ ) 
has singularities for the solutions of (^.6|) and, moreover, there is a special point c = — r^^ at 



which the function 0(c) becomes divergent. As was shown in sect. 5, for the solutions of (|5.6| ) 



the string susceptibility exponent is equal to "fstr = 1/2. Being combined together, equations 
(|5.6| ) and ( [4.3|) define the curve of criticality at the {a,g) plane with the exponent jstr 



1 

2' 

However, in despite of the model with polynomial potential considered in the previous sections 
this phase cannot be identified as the phase of branched polymers because it has the following 
"anomalous" properties. Firstly, the cosmological constant acr{g) decreases in this phase when 
the touching constant increases. Recall, that at the phase of branched polymers, considered 
in sect. 7, acr{g) has an opposite behaviour shown in Fig. 1(a). Secondly, at the end point of 
this curve with Ocr = —2 and g = A/r]'^ the parameter Ccr calculated from (|5.6|) turns out to be 
equal to Ccr = —ff giving rise to the singularities of f{ccr) and 0(ccr). As a consequence, the 
derivative ^ defined in ( |5.8|) acquires additional singularities leading to the increasing of the 



critical exponent at that point to 'jstr = 2. This number lies outside the allowed region of values 
7str < 1/2 one may get in any model of noninteracting planar surfaces with positive weights 
TB| . Moreover, the effective touching constant scales near this point as (7 ~ 1/ y^(a — ctcr) and 
tends to infinity as a ^ a^r- Thirdly, the phase diagram of the model has two different curves 
corresponding to 7^4^ = and 7^4^ = 2. It turns out that despite of the analogous diagram 
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for the perturbed matrix model with polynomial potential shown in Fig. 1 these curves do not 
intersect each other. All these properties differs from what one expects to get by considering 
the perturbed Penner model as discretization of a sum over surfaces for Polyakov's bosonic 
string. This means that the Penner model originally proposed to capture Euler characteristics 
of the moduli space of Riemann surfaces fails to describe random surfaces after perturbation 
by higher order curvature term. 

Note added 

After the completion of this paper 1 was informed by L. Alvarez-Gaumc that he got analogous 
results applying large reduced models to the study of string with c > 1. 
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